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The presence of a space charge region induces an internal electric field within the charged region 
that, in a ferroelectric material, would rotate the polarisations to align with the field. The strength 
of the indnced field would therefore determine the domain patterns and polarisation switching 
properties of the material. Using a phase-field model, we investigate the effect of charged layers in 
fully and partially depleted BiFeOs thin films in the rhombohedral phase. While the domain pattern 
in a charge-free BiFeOs film consists of only two polarisation variants, we observe complex patterns 
with four coexisting variants that form within the charged layers at sufficiently high induced fields. 
These variants form a head-to-head configuration with an interface that is either wavy or planar 
depending on the internal field strength, which is determined by the charge density as well as the 
thickness of the charged layer. For depletion layers with sufficient thickness, there exists a range of 
charge density values for which the interface is wavy, while at high densities the interface becomes 
planar. We find that films with wavy interfaces exhibit enhanced susceptibilities with reduced 
hystereses compared to the charge-free film. The results of our work suggest that introducing space 
charge regions by careful selection of dopant density and electrode materials can engineer domain 
patterns that yield a higher response with a smaller hysteresis. 


I. INTRODUCTION 

When a metal electrode is placed in contact 
with a semiconductor material, there may be a 
flow of charge carriers from the semiconductor 
to the metal, depending on the work function of 
the metal [1]. An accumulation of charge carri¬ 
ers results at the electrode/semiconductor inter¬ 
face, inducing an electrostatic potential that de¬ 
pletes these carriers in adjacent regions [2, 3]. This 
space charge region, or depletion layer, is typically 
formed near interfaces such as grain boundaries 
and interfaces between heterophases, as well as 
near surfaces. In ferroelectric materials which are 
regarded as wide band gap semiconductors, the 
presence of such regions may affect the properties 
of the material, in many cases detrimental for ap¬ 
plications such as non-volatile memory devices and 
capacitors. It lowers the ferroelectric phase tran¬ 
sition temperature and suppresses ferroelectricity 
in the space charge regions [4]. Furthermore, the 
domain pattern that forms under the influence of 
space charge can substantially differ from those 
without [5-8]. Consequently when such a sample 
is switched, the hysteresis loop constricts with a 
lowering of the coercive field, and depending on 
the extent of the depletion layer, the remnant po¬ 
larisation as well [2, 8-10]. If the distribution of 
charges in the sample is asymmetric, the result¬ 
ing internal electric field leads to a preferred po¬ 
larisation direction, thereby shifting the hysteresis 
loop along the electric field axis; this phenomena 
is called imprint [8, 11]. 

Using a scanning probe microscope (SPM) tip 
to apply an electric field over a small region in 
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a (100) BiFeOs film) Vasudevan et al. observed 
the formation of domains with polarisations that 
form a closed loop (i.e. head-to-tail configuration 
of domains) at low applied fields, and domains 
with polarisations directed towards a central point 
(head-to-head configuration) at very high fields. 
[12]. Similar patterns may also be formed within 
space charge regions in ferroelectric films due to 
the internal electric field induced by the charges. 
Head-to-tail and 180° domains, whose walls are 
charge-neutral, are the more stable configurations 
in a charge-free material, whereas charged walls 
between head-to-head or tail-to-tail domains may 
stabilise in the presence of space charge [13-16]. 
Usually in theoretical investigations concerning 
the latter type of domains, the models have flat 
interfaces between the domains [15-18]. Other 
studies, both experimental [19-21] and theoreti¬ 
cal [22] , report the formation of zigzagged domain 
walls. It has been suggested that the formation of 
such interfaces allows the electrostatic energy to be 
reduced, as locally at the domain walls the polari¬ 
sations are rotated into head-to-tail configurations 
[19]. Misirlioglu et al. find a critical thickness of 
the space charge layer above which this type of 
domain wall develops in their simulated BaTiOa 
films, and below which the film remains as a sin¬ 
gle domain [22]. 

It can therefore be deduced that an approach 
to engineer complex domain structures is to ma¬ 
nipulate the space charge regions in ferroelectric 
materials. 

The thickness of space charge layers in thin films 
is determined by the built-in voltage across the 
layer. For a given dopant concentration, the built- 
in voltage is a property that depends on the ma¬ 
terial of the electrode. The layer thickness, which 
determines whether the film is fully or partially 
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depleted, is thus dependent on the selection of 
dopant concentration and the electrode in experi¬ 
mental set-ups. When simulating a charged layer 
within a film, the layer thickness and dopant con¬ 
centration may be chosen as variables instead; this 
would be equivalent to specifying the built-in volt¬ 
age across the layer. 

In investigating the properties of ferroelectric 
materials, the phase-field method has been shown 
to be a powerful technique and has been widely 
used [23]. A real-space approach with this tech¬ 
nique takes into account long-range elastic and 
electrostatic interactions in the material of study, 
and is convenient for nanoscale systems in which 
free surfaces are important since it is straightfor¬ 
ward to apply the boundary conditions [24]. This 
method has been applied in simulating the domain 
patterns in free-standing nanostructures of tetrag¬ 
onal lead zirconate-titanate [25]. 

Most of the studies of space charge effects in 
ferroelectric materials have been devoted to mate¬ 
rials such as BaTiOa [6, 16, 17, 21, 22], PbTiOs 
[8, 13, 15] and PbZri_xTi ,,03 (PZT) for x > 0.47 
[10, 11], all of which adopt the tetragonal ferro¬ 
electric phase. These systems have 90° and 180° 
domain walls, with only the former producing a 
ferroelastic configuration. In these works, the ex¬ 
ternal electric field is usually applied parallel to 
the polar direction. 

Few have investigated the effect of space 
charge on rhombohedral ferroelectric materials 
[12, 19]. Unlike tetragonal ferroelectric materi¬ 
als, the rhombohedral phase has a larger num¬ 
ber of ferroelectric variants with two possible fer¬ 
roelastic domain patterns of 71° and 109° domain 
walls. Thus it can have more complex domain pat¬ 
terns [26] , whose piezoelectric response may be en¬ 
hanced especially when applying an electric field 
that is not along one of the polar directions [17, 27- 
32]. 

A rhombohedral ferroelectric material, BiFeOa, 
has been attracting much interest recently be¬ 
cause of its multiferroic [33, 34] and photoelectric 
[35, 36] properties. It has potential applications in 
memory storage, solar energy storage, sensors and 
telecommunications devices [34, 37-39]. 

Therefore in this work, we select BiFeOa with 
polarisations along the (111) directions (with ref¬ 
erence to the pseudocubic lattice) as our mate¬ 
rial of study and examine the domain structures 
that form in [001] thin films under the influence 
of space charge. We also investigate the effective 
polarisation-electric field behaviour when applying 
an electric field in the [001] direction. 

Assuming static charges on the timescale of po¬ 
larisation switching [6] for the case where the space 
charge region is generated due to the presence of an 
electrode/film interface, what we find is a wealth 
of domain patterns that varies with the strength 
of the electric field induced by the charges. 

Calculating the response in the direction par¬ 


allel to the applied field, we note an optimum 
charge density with minimum charged layer thick¬ 
ness that produces a domain pattern with a wavy 
interface whose mobility is higher than the others. 
Our findings would be useful for applications that 
require materials with high piezoelectric response 
and small hystereses. 


II. COMPUTATIONAL DETAILS 


The dynamics of the polarisation order parame¬ 
ter components Px, Py and Pz evolving with time 
t is given by the time-dependent Ginzburg-Landau 
equation as 


dt 


= -r 



( 1 ) 


where T is a kinetic coefficient related to the do¬ 
main wall mobility, F is the free energy of the 
ferroelectric material and Ei is the component i of 
the electric field E in the sample. F consists of the 
Landau free energy Fi , the gradient free energy as¬ 
sociated with the domain walls Eg, the elastic free 
energy Eg and the electrostatic free energy E^i: 


F — Fl + Fg + Ff. + Fei (2) 

We express Fi in terms of Px, Py, and Pz up to 
the fourth order: 


Fi = 


Jdr ai {Pi 


F Py F Pz) F an {P^ + P^ + P^) 


+ an [PlPl + PlPl + Pi Pi) 


( 3 ) 


where ai, an and ai 2 are material parameters. 
Fg can be expressed as: 

Fg=J dr [{VPxf + {VPyf + {VPzf 

( 4 ) 

where AT is a gradient coefficient, while the expres¬ 
sion for Fe is as below: 


Ff, = dr 
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Cii, Ci 2 and C 44 are material constants. Sij are 
the strain tensors having the usual definition in 
terms of displacements Ui 


1 / dui duj \ 

2 dxi) 


( 6 ) 


and the electrostrictive strain tensors e?- have the 
following expressions: 

el=QiiPf + Qi2{P] + Pl) ( 7 ) 
and =Qi 4 PiPj for i ^ j ( 8 ) 

with Qii, Qi 2 and Q 44 being material-specific elec¬ 
trostrictive constants. 

We determine the electric field E in the sample 
from the electrostatic potential (/) using the rela¬ 
tionship 


E = -V0 (9) 

as well as Gauss’s law relating the electric displace¬ 
ment D, E and polarisation P: 

V-D = V-(eoehrE -t- P) = -eoebrV^(/)-|-V-P = p (r) 

( 10 ) 

Cq is the permittivity of free space, e^r = 10 is a 
background permittivity [40] and p (r) we have set 
as 


Piz) 


po, L - z <w 
0 , otherwise 


( 11 ) 


where L is the thickness of the film and w is the 
thickness of the space charge layer. We assume 
Po — qNjj where q is the charge of a carrier and 
Nu is the defect (or carrier) density. In our simu¬ 
lations, we arbitrarily select q = —e where e is the 
elementary charge. We now have Fei, defined as 

Pel = J dr • E,) ( 12 ) 


The dissipative force balance equations give the 
displacement field dynamics: 


d'^Ui 





duj 

dt 


(13) 


where p is the density and p is a viscous damping 
term. The latter functions to drive the system to 
mechanical equilibrium such that = 0. The 
stress tensors aij are obtained as 


5F 

5Sij 


( 14 ) 


We consider a (001) epitaxial BiFeOa film on 
top of a substrate, with an electrode deposited at 
the top surface of the film, and simulate the dy¬ 
namics of the polarisations by numerically solving 


TABLE I. The material constants for rhombohedral 
BiFeOs taken from reference 41 . T is the temperature 
in Kelvin. 

«! 

4 . 9 (T- 1103 ) X 10 ^ V m/C 

On 

6.5 X 10 ® V mVC® 

012 

1.0 X 10 ® V mVC® 

Cii 

3.00 X 10 “ N/m 2 

C12 

1.62 X 10 “ N^/m^ 

C44 

0.691 X 10 “ NVm 2 

Qii 

0.032 m^/C^ 

Qi2 

- 0.016 m^/C^ 

Q44 

0.02 m^/C^ 


the model described in equations 1 to 14 using fi¬ 
nite differences. Initialising P^, Py and Pz with 
small, random fluctuations about zero to simu¬ 
late the paraelectric initial conditions, we exam¬ 
ine the ferroelectric domain patterns that form in 
fully and partially depleted films. We employ a 
128nm x 128nm x 32nm simulation cell with peri¬ 
odic boundary conditions for P, u and (p in the lat¬ 
eral directions. The top surface is stress-free with 
aij ■ n = 0, where n is the surface normal, and the 
bottom surface is clamped {ux = Uy = Uz = 0 ). 
The film is also laterally clamped as we have pe¬ 
riodic boundary conditions on Ux and Uy later¬ 
ally. Assuming that the top and bottom interface 
charges compensate exactly the respective bound 
interface charges, we impose = 0 at z = 0 and 
z = L = 32 nm. Rescaling all lengths as r = r*(5, 
we choose d = 1 nm as our grid spacing. We rescale 
the polarisation as P = EqP* where Eg = 0-9 
Cm“^ is the spontaneous polarisation of BiFeOg, 
the strains as = eoe*j where eg = QiiEg; and 
the time as t = t^t* where tg = (E|ai|)“^. All 
other coefficients are consequently rescaled in the 
same way as reported in reference 25, and we use 
their values for rescaled coefficients of E, E, p and 
p. The material specific parameters are taken from 
reference 41 and listed in Table I. We model films 
with values of Nd between 10 ^® m“^ to 10 ^^ m“^, 
which are within the range of those for heavily 
doped BiFeOg films [42, 43]. The charges are uni¬ 
formly distributed within space charge layers par¬ 
allel and adjacent to the film/electrode interface 
with thicknesses w of 0, 5 , 10, 15 and 32 nm, with 
the latter thickness corresponding to the fully de¬ 
pleted case. 


III. RESULTS AND DISCUSSIONS 

A. The charge-free BiFeOa thin film 

Bulk BiFeOg in the rhombohedral phase has 
eight possible polarisation variants along the ( 111 ) 
directions, which are commonly denoted as rf = 
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FIG. 1. (a) The domain pattern of BiFeOa thin film without presence of charge, with the direction of polarisation 
indicated for y — 0 plane. The colour within the film distinguishes between the domains. This domain pattern 
also corresponds to the point labelled (ii) in (c). (b) The domain pattern of BiFeOs thin film without presence 
of charge corresponding to the remnant state labelled (i) in (c). (c) The simulated Pz — Eapp hysteresis loop for 
the film. 


[in], r+ = [Ill], r+ = [III], r+ = [ill], = 
[111], r 2 = [111], ^3 = [111], r 4 = [ill] (with ref¬ 
erence to the pseudocubic lattice). Without any 
presence of charges in the film, our BiFeOa thin 
film evolves from the paraelectric state to form 
71° twinned domains with domain walls at an an¬ 
gle greater than 45° to the film/substrate inter¬ 
face, consistent with results reported by Zhang et 
al. using nearly similar material constants in their 
phase-field model [44]. Figure 1(a) shows 
domains that form in some of our simulations when 
steady state is reached. With different random 
seeds to initialise the polarisations in the paraelec¬ 
tric state, we have observed other twinned domain 
formations with {101}-type walls as suggested by 
Streiffer et al. for films with rhombohedral lattices 
[45] such as r'^/rf, and r^ , to name a 

few. They also noted that all domains in a film 
must have the same Pz direction to maintain a 
charge-neutral film. In all of our simulated films, 
the domain walls deviate slightly from the ideal 
{101} orientation; Zhang et al. have shown that 
this deviation is due to the competition between 
the domain wall energy that minimises when the 
domain wall is perpendicular to the film/substrate 
interface, and the elastic energy that favours the 
ideal orientation at 45° to the interface [44]. 

We then apply a varying potential difference 
across the thickness of the film with /r^ do¬ 
mains to simulate the ferroelectric switching of the 
charge-free film, by varying the external electric 


field from 0 kV/cm to 3000 kV/cm, —3000 kV/cm 
and back to 0 kV/cm in t* = 10^ timesteps. Com¬ 
puting the average Pz as a function of the applied 
field Eapp we show the resulting curve in Figure 
1(c). The r]”/r{i domains are observed for the 
remnant state with Pz < 0, corresponding to the 
point labelled as (ii) in Figure 1(c). These domains 
switch fully to /r^ domains upon forward bias, 
and back to /r~^ domains with reverse bias. The 
/r}" domain pattern, illustrated in Figure 1(b), 
is stable at the remnant state with Pz > 0 which 
is labelled as (i) in Figure 1(c). The P^ and Py 
components of the polarisation in the film remain 
unchanged during the switching as we expect, with 
only the Pz component switching, since we are ap¬ 
plying the electric field in the [001] direction. As 
we switch between the two remnant states, there is 
no long-range migration of the domain wall; only 
the inclination of the wall is altered. 


B. Domain formation in fully and partially 
depleted films 

When we include uniformly distributed static 
charges throughout the film to model a fully de¬ 
pleted film, we find that instead of having only 
twinned r^/r^ domains in the film, the top of 
the film consists of a set of [Px,Py,—Pz] and 
[Px,—Py,—Pz] twinned domains and the bottom 
half, [Px,Py,Pz] and [Px,-Py,Pz] domains. Fig- 
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{si)No = 3 X 10^® m-3 (h)ND = 6 x 10^® (c}Nn = 1-1 x 10^’’ 

FIG. 2. Domain patterns in fully depleted films with No values from 3 x 10^® m“® to 1.1 x 10^^ m“®. Lighter, 
brighter shades represent domains with more positive Pz (i.e. r+) and darker shades, those with more negative 
Pz (i.e. r“). 



FIG. 3. The profile of the internal field Ez (in the 
paraelectric phase) induced by uniformly distributed 
charges in fully depleted films with Nn values ranging 
from 3x 10^® m“® to l.lx 10^^ m“® across the thickness 
of the films. 


ure 2 illustrates these domain configurations form¬ 
ing in fully depleted films for different values of 
N]:,, clearly demonstrating that the geometry of 
the domain wall between these sets of twins de¬ 
pends on the magnitude of the charge density. The 
film with a lower value of Njj = 3 x 10^® m“^ 
has a wavy domain wall in the lateral directions 
of the film separating the sets of twinned domains 
that meet in a head-to-head configuration (Figure 
2(a)). For high values of Njj, this wall has a planar 
geometry (Figures 2(b) and 2(c)). 

To understand the formation of the domain con¬ 
figurations described above, we examine the inter¬ 
nal electric field induced by the charges. A uni¬ 
form charge density produces a radial electric field 
inside a charged body. For the thin film geome¬ 
try considered in our work, the field varies only 
in the z direction. This internal field may influ¬ 
ence the domain pattern that forms when the sam¬ 
ple is cooled from the paraelectric phase to below 
the ferroelectric transition temperature. Since the 
field is already present in the paraelectric phase, 
we perform additional simulations in the paraelec¬ 


tric phase for the same values oi Njj as before to 
determine the variation of the internal electric field 
within the films. As shown in Figure 3, the internal 
electric field generated due to the uniform distri¬ 
bution of charges in a film is equal in magnitude 
but opposite in direction in the top and bottom 
halves of the film. Since we have introduced neg¬ 
ative charges in the film, the electric field in the 
bottom half is positive while that in the top half is 
negative. The magnitude of the field is maximum 
at the top and bottom surfaces and decreases to 
zero towards the middle of the film. 

The electric field will induce a rotation of the 
polarisations to align with the field locally if the 
magnitude of the field is large enough. For films 
with sufficiently high values of Nd to generate 
large enough internal fields, the positive field in 
the bottom half of the charged region stabilises 
domains with > 0 and the negative field in the 
top half, domains with Pz < 0. In rhombohedral 
ferroelectrics, this means that only rf, r^, and 
rf may form in the bottom half and r(~, 
and rj in the top half. In our simulations, two of 
the four possible domains in each half are formed, 
where the domains in one half must each be twins 
of those in the other half through the (001) twin 
plane, giving rise to the four-domain pattern ob¬ 
served in Figure 2. Thus it is the inhomogeneity in 
the internal electric field that is responsible for the 
formation of multiple domains within the charged 
region. 

Comparing the variation of the internal field 
for films with different values of Nd (Figure 3), 
we observe that the magnitude of the internal 
field within the charged region increases as Al^ 
increases. For low values of the small magni¬ 
tude of the internal field is insufficient to produce 
a pattern of domains different from that in the 
charge-free film. For depleted films with interme¬ 
diate values of Nu, the magnitude of the internal 
field which is larger in the near-surface regions, is 
sufficient to induce the formation of domains with 
polarisations that align with the field in these re¬ 
gions. Towards the middle of the film, however, 
the magnitude of the field is too small to force the 
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FIG. 4. Polarisation vector maps of Pz — Py in an 
X — z plane of fully depleted films with No values of 
(a) 3 X 10^® ni~^, (b) 6 x 10^® m“® and (c) 1.1 x 10^^ 
m“®. Domains rj (top left of film), rj (bottom left 
of film), r)” (top right of film) and rj (bottom right 
of film) are indicated. These maps correspond to the 
domains illustrated in Figure 2. 


polarisations to rotate with the field. In this case, 
the polarisations rotate into a head-to-tail configu¬ 
ration locally near the domain wall, accommodat¬ 
ing the constraint imposed by the induced domains 
in the near-surface regions and at the same time 
minimising the electrostatic energy. Consequently 
the interface takes on a wavy geometry, as illus¬ 
trated in Figure 4(a) showing the — Py vectors 
in an a; — z plane for the fully depleted film with 
Nd=3x 10^® m-3. 

When the value of Njj is high, the magnitude 
of the electric field increases more rapidly from 
the middle of the film. The field is therefore large 
enough in the middle of the film to induce a head- 
to-head domain configuration with a planar in¬ 
terface (Figures 4(b) and (c)). It also results in 
a rotation of the dipoles from pure rhombohe- 
dral orientations towards the [001] and [OOl] di¬ 
rections, such that for the fully depleted film with 
No = 1.1 X 10^^ m“^, the domains approach the 
tetragonal phase. The increasing P^ component 
in the near-surface regions of films with increasing 
No can be observed from the Pz — Py vector maps 
in Figure 4, and we also represent the degree of 
the polarisation rotation with the changing shades 
within the domains in Figure 2. 

In partially depleted films, the presence of a 
charged layer shows a similar effect. For films with 
uniformly distributed charges of No = 3 x 10^® 
m“® and No = 6 x 10^® m“® within charged layers 
of width w = 5 nm at the top surface of the films, 
the domain patterns resemble that of a charge- 
free film with twin /r^ domains as shown in 


Figures 5(a) and (b), except for a slight decrease 
in Pz near the top of the layer (indicated by the 
different shades in this region). When we increase 
No to 1.1 X 10^^ m“® for the same value of w, we 
observe the formation of the twinned set of 

domains in the top half of the charged layer and 
the r'^/r^ set in the bottom half of the charged 
layer, with the latter set extending into the charge- 
free region (Figure 5(c)). 

We obtain a similar trend for films with w = 10 
nm, but with the appearance of the double set of 
twinned domains from lower values of No- Small 
regions of /r^ domains already form at the do¬ 
main walls in the film with No = 3 x 10^® m“^, 
represented as the dark-shaded regions in Figure 
5(d). Increasing w to 15 nm with No = 3 x 10^® 
m“® leads to the appearance of such domains in 
a larger number (Figure 5(g)), such that the do¬ 
main pattern begins to resemble that of the fully 
depleted film having the same No with the wavy 
domain interface (Figure 2(a)). The /r^ set of 
domains from the lower half of the film extends up 
to the top surface of the film in some regions, as 
w = 15 nm is not sufficiently large to allow the for¬ 
mation of a continuous wavy interface. For films 
with this value of w, and No values of 6x 10^® m“® 
and 1.1 X 10^^ m“®, the two sets of twinned do¬ 
mains form with planar interfaces, with rotation of 
the polarisations approaching the tetragonal phase 
as before even within the charge-free regions. 

Whether the set of r]" /r^ domains forms within 
the charged layer of a partially depleted film de¬ 
pends on the strength of the internal electric field, 
determined by the values of w and No- The field 
within the charged region, whose profile is shown 
in Figure 6 for films with No = 6 x 10^® m“®, 
does not induce the formation of these domains 
in the top half of the charged layer if the magni¬ 
tude of the field is too small. As the magnitude 
increases with increasing w and No^ the do¬ 

mains are observed near the top of charged layers 
with higher values of these parameters. Within 
the much larger charge-free region, the internal 
field is non-zero and positive, increasing in mag¬ 
nitude with increasing No and w- Although the 
magnitude in this region is much smaller than the 
maximum magnitude attained at the top of the 
charged layer, it is sufficient to bias the formation 
of the /r^ domains within the charge-free re¬ 
gion, with the domains approaching tetragonality 
with increasing magnitude of the internal field. 

The appearance of the domain pattern with the 
wavy domain interface for ic > 15 nm is consis¬ 
tent with the finding by Misirlioglu et al- From 
their 2-dimensional Ginzburg-Landau-Devonshire 
model of depleted tetragonal BaTiOa films, they 
predict a critical thickness above which head-to- 
head domains with a zigzagged interface develop 
[22]. We also observe a critical value of w associ¬ 
ated with the strength of the induced electric field 
for the appearance of the domain pattern with pla- 
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FIG. 5. Domain patterns in partially depleted films with depletion layer thicknesses from 5 nm to 15 nm, and 
No values from 3 x 10^® m“® to 1.1 x 10^^ m“®. 


X 10* 



FIG. 6. The Ez profile across the thickness of partially 
depleted dims with No = 6 x 10^® m“®, for values of 
w from 5 nm to 15 nm. 


nar horizontal interface that we find in films with 
high Nd. For films with N^ = 6 x 10^® m“®, this 
configuration appears when ic > 10 nm (Figure 
5(e)), while films with No = 1 x 10^^ m“® need 
only w > 5 nm for the same configuration to form 
(Figure 5(c)). 


There is in addition a critical value of No be¬ 
yond which the wavy domain wall transitions to 
a planar interface. We have observed this pro¬ 
gression for the fully depleted films, as well as for 
the partially depleted films with w = 10 nm (not 
shown here) and 15 nm (Figures 5(g-i)). For a film 
with a thinner depletion layer of w = 5 nm, the 
charged layer is too thin to allow the formation 
of a wavy wall. Thus the twinned domain pattern 
resembling the charge-free film directly transitions 
to the four-domain pattern with planar interface 
when No > 1.1 X lO^^m"®. 

Various experimental studies have reported ob¬ 
servations of domains in a head-to-head or tail- 
to-tail configuration with a wavy interface, one of 
which have recently been observed in a tetragonal 
PZT film [20]. The authors in that work attributed 
the domain shapes to the anisotropy of domain 
growth, and suggest opposite internal fields near 
the top and bottom surfaces of the film due to 
high oxygen vacancy concentration and electronic 
band bending at the film/substrate interface re¬ 
spectively. These are thought to induce domains 
with polarisations pointing down towards the mid¬ 
dle of the film (Pz < 0) near the top of the film, 
and domains with upwards polarisations {P^ > 0) 
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at the bottom of the film. Our results show that 
domains in a head-to-head configuration may be 
generated by a single negatively charged layer if 
the magnitude of the induced internal field is suf¬ 
ficiently high, with a residual field even within the 
charge-free regions. Furthermore, the geometry of 
the interface is dictated by the strength of the 
field. Having two oppositely charged layer in a 
film may produce the same domain configuration 
provided the induced field within the other posi¬ 
tively charged layer is small enough in magnitude. 
If the magnitude is sufficiently high, additional do¬ 
mains in tail-to-tail configuration would result in 
the positively charged layer. 

We observe that in both partially and fully de¬ 
pleted films that we have considered in our work, 
the induced internal fields are very large in mag¬ 
nitude. In principle, such large fields in the ferro¬ 
electric state may lead to charge transport which 
can modify the field and thus the domain pat¬ 
terns, particularly in films with large Nij. How¬ 
ever, we expect the charge migration to be a slow 
process. For example, in a study of the ageing 
phenomenon in ferroelectrics, it has been shown 
that charge migration induced by a depolarisation 
field leads to a build-up of space charge over time 
scales ranging from 10® seconds to 10^ seconds [46]. 
This is much slower than the time scales associ¬ 
ated with domain reorientation, which is on the 
order of 10“® seconds. Thus we expect findings 
based on a model with static space charge to be 
valid for short times and fast switching frequen¬ 
cies. The fact that the resulting domain patterns, 
such as those with wavy interfaces, have indeed 
been observed in experiments for various systems 
[19-21] gives us confidence that our observations 
will still hold in the long time limit, as the internal 
field may not be completely eliminated by mobile 
charges. 

C. Polarisation switching and response of 
the films 

When an external electric field is applied to a 
film with depleted layers, the internal field pro¬ 
duced by the presence of charges may oppose the 
effect of the externally applied field. Within the 
thin films that we have considered in this work, 
we have shown that the internal field along the 
^-direction is positive in the lower regions of the 
charged layers but negative in the upper regions 
in the absence of external bias. Therefore one of 
these regions will always oppose an applied field 
parallel to the internal field. Consequently, the do¬ 
mains will be pinned, requiring prohibitively large 
applied fields to completely switch to domains hav¬ 
ing the same polarity in z. 

Figure 7 demonstrates these effects for the fully 
depleted film with a wavy domain wall {Nd = 
3 X 10^® m“®) under the same switching condi¬ 


tions that we have applied for the charge-free film 
in Figure I. We also illustrate the domain patterns 
at maximum magnitude and zero applied fields in 
the same figure. The magnitude of the average Pz 
increases with the magnitude of the applied field, 
and the wavy domain wall moves towards the top 
or bottom surface of the film. However, not all 
regions of the domains switch to align with the ex¬ 
ternal field as evident from the presence of darker 
shaded r~ domains remaining at the top surface 
of the film in Figure 7(a) and lighter shaded r+ 
domains at the bottom surface of the film in Fig¬ 
ure 7(c). The internal field induced by the charges 
at the top of the film opposes the external field 
under forward bias, and that at the bottom of the 
film when the field is reverse-biased, leading to in¬ 
complete switching of domains. The magnitude 
of average Pz at maximum magnitude of applied 
field is thus less than the saturation polarisation 
of a charge-free film. 

When the external field is removed, the internal 
field restores the four domains almost to the orig¬ 
inal state, where the average Pz is close to zero. 
Therefore the remnant Pz is close to zero and the 
magnitude of an applied field required to reduce 
the magnitude of Pz to zero, i.e. the coercive field, 
is also close to zero; correspondingly, the hysteresis 
loop narrows. This constriction of the hysteresis 
loop, as well as the decrease in the remnant polar¬ 
isation, has been observed in earlier works [2, 10]. 

The hysteresis loops for the fully depleted films 
become more constricted with increasing values of 
Nu, as shown in Figure 8(a). This is consistent 
with the trend observed in other works [8, 9]. The 
increasing magnitude of the internal field in the 
films increasingly opposes the polarisation switch¬ 
ing, such that for films with No > 6 x 10^® m“® 
there is no switching of domains but only a slight 
enhancement in the magnitude of average Pz with 
Eapp- With very little movement of the domain 
walls in the z direction, the hysteresis completely 
disappears. The curve becomes nearly horizontal 
and linear, reminiscent of that for a paraelectric 
material. 

We also observe a similar trend for partially de¬ 
pleted films, as shown in Figure 8(b) for films with 
w = 5 nm. In particular, we find that for the par¬ 
tially depleted film with No = 1.1 x 10^^ m“® 
which forms the four-domain pattern in the ab¬ 
sence of external field (Figure 5(c)), the rf/rj do¬ 
mains at the top of the charged layer do not switch 
to r^/r 2 domains under forward bias (Figure 
9(a)), although complete switching of the r^/r^ 
domains to jr^ domains occurs with reverse 
bias (Figure 9(b)). The field induced within the 
charged layer in this film follows the same trend 
as those shown in Figure 6: the internal field at the 
top of the charged layer, where Ez < 0, is much 
larger in magnitude than elsewhere in the film. It 
is sufhciently large to impede the switching of the 
r]"/r^ domains, whereas the induced field within 
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FIG. 7. Simulated hysteresis loop (left) of fully depleted film with No = 3 x 10^®m and the domains that 
form (right two columns) at points (a)-(d) during switching. 



(a) 



FIG. 8. Simulated hysteresis loops of (a) fully and (b) 
partially depleted films with w = 5 nm, for No values 
between 1 x 10^®m“®, to 1.1 x 10^^m“®. The curve for 
the charge-free film is included for comparison in both 
graphs. 
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FIG. 9. Polarisation vector maps of Pz—Py in an a; — 2 : 
plane of the partially depleted film with w = 5 nm 
and No = 1.1 x 10^^ under an external electric field 
Eapp at (a) forward bias and (b) reverse bias. The 
r J (top left of film) and r)” (top right of film) domains 
indicated in (a) are pinned by the internal field induced 
within the charged layer under forward bias, whereas 
the rj (bottom left of him) and r'^ (bottom right of 
him) domains switch completely to r J and rj” domains 
respectively with reverse bias. 


the r'^/r^ domains is too small to pin the polar¬ 
isations. This incomplete switching of domains, 
which we also observe for other partially depleted 
films with the head-to-head domain configuration, 
leads to the observation of non-switchable regions 
in ferroelectrics [20]. 

There is in addition a shift of the hysteresis loop 
along the Eap-p axis that is not present in that of a 
fully depleted film, commonly referred to as im¬ 
print. Due to the asymmetrical distribution of 
the charges, the upper region of the charged layer 
where the internal field E^ < 0, is much smaller 
than elsewhere in the film where E^ > 0 (Figure 
6). In the absence of an external field, the av¬ 
erage value of the internal field is positive in the 
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FIG. 10. x/Xo of partially and fully depleted films as 
a function of No- 


z direction despite the presence of the local neg¬ 
ative field of much larger magnitude, and hence 
domains with > 0 will be the preferred polar¬ 
isation state. Consequently, a larger electric field 
is required to switch the domains such that an av¬ 
erage Pj; < 0 is produced in the film, compared to 
that for the reverse. The resulting hysteresis loop 
is thus left-shifted and becomes further left-shifted 
with increasing strength of the internal field, as 
can be observed for the film with w = 5 nm in 
Figure 8 (b). 

Having demonstrated the influence of charges 
on domain patterns and polarisation switching be¬ 
haviour, we now consider the implications on the 
dielectric response. From the P^ — Eapp curve, we 
calculate 3 ^^ at zero applied field for P^ > 0 , 
which gives an indication of the susceptibility of 
the film in z, y = Figure 10 shows this 

quantity relative to tliat of the charge-free film, 
x/Xoj as a function of Njj for partially and fully 
depleted films. 

We find that an initial increase in No results 
in increasing enhancement in the susceptibility of 
the films. The magnitude of the internal field 
at top/bottom surfaces of a charged layer, when 
large enough, rotates and pins the polarisations 
from P^ > 0 to Pz < 0 in the top region of the 
charged layer, thus leading to the appearance of 
head-to-head domains that are stable even at re- 
manence; in comparison, a charge-free film at re- 
manence would have a single polarisation direc¬ 
tion in z. Figure 11 illustrates the increasing vol¬ 
ume fraction of r)~ /domains at remanence with 
Njj. As 180° switching of the local polarisations 
Pz in films with head-to-head domains produces 
a larger polarisation response compared to that 
of domains with a single polarisation direction in 
z, the susceptibility of the film therefore enhances 
with Njj. Additionally we observe that films with 
higher values of w, and thus having larger volume 
fractions of domains, show greater enhance¬ 

ment in susceptibility. The response of the fully 


depleted film exhibits an increment of almost four¬ 
fold when Nd = 2.3 x 10^® m“^, relative to that of 
the charge-free film. The partially depleted films 
with It; = 15 nm also show significantly enhanced 
response, though lower than that of the fully de¬ 
pleted film, of almost twofold when No = 4 x 10^® 
m“^. These films have in common wavy inter¬ 
faces between domains in head-to-head configu¬ 
ration, with hysteresis loops that narrow without 
complete loss of hysteresis and without significant 
shifts along the Eapp axis. 

Below these No values, the internal field that 
acts to pin the domains in the head-to-head con¬ 
figuration is relatively weak in the vicinity of the 
domain wall, allowing switching of the polarisa¬ 
tion with externally applied fields. However, there 
is a competing effect also as a result of increas¬ 
ing internal field strength with No that inhibits 
the polarisation switching. The movement of the 
domain walls becomes more limited and the sus¬ 
ceptibility then decreases with further increase in 
No, eventually falling to a value that is lower than 
that of the charge-free film, as can be observed in 
Figure 10. 

Furthermore, the maximum in the susceptibil¬ 
ity occurs at lower No values with increasing w, 
as the magnitude of the internal field that inhibits 
domain switching also increases with w. For films 
with suppressed susceptibilities, we have observed 
that the interface becomes planar and the hys¬ 
teresis loop becomes almost linear or further left- 
shifted. 

Thus there exists a range of values for No for 
a given value of w that enhances the susceptibil¬ 
ity of a film containing a depletion region, beyond 
which is detrimental to the susceptibility. This 
indicates that the presence of space charge in fer¬ 
roelectric materials may be useful for applications 
that require high susceptibilities as well as small 
hystereses. The enhancement due to the presence 
of charges is in addition to the already known ef¬ 
fect of poling along a non-polar direction in a fer¬ 
roelectric material. 


IV. CONCLUSIONS 

We have investigated the influence of space 
charge on domain patterns in rhombohedral ferro¬ 
electric thin films, taking BiFeOs as an example. 
Phase field simulations of domain patterns and po¬ 
larisation switching have been performed for differ¬ 
ent charge densities and charged layer widths. We 
show that the introduction of a charged layer cre¬ 
ates an internal electric field which may lead to the 
formation of double twinned domains in a head-to- 
head configuration. Depending on the magnitude 
of the internal field which we vary with charge den¬ 
sity and layer width, we find a wavy or planar in¬ 
terface between the head-to-head domains. The 
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FIG. 11. Domain patterns at remanence {Pz > 0) of fully depleted films with No values from 1.1 x 10^® m ® to 
2.3 X 10^® m“^, showing the increasing amount of jr^ domains forming with No- 


former forms in films with intermediate values of 
Njj above a critical value of w, while the latter is 
observed when No is large. 

We have also examined the polarisation switch¬ 
ing behaviour of the fully depleted as well as the 
partially depleted films with the applied field along 
the [001] direction. The widely observed reduction 
in saturation and remnant polarisation, reduction 
in coercive field, constriction of the hysteresis loop 
and, for partially depleted films, its translation 
along the applied electric held axis, is explained in 
terms of the internal held that pins or imprints the 
domain structures as they evolve with the applied 
switching held. For hlms in which the charged 
layer induces a sufficiently large held, the pinning 
of the domains also results in regions with non- 
switchable polarisations. We hnd that the suscep¬ 


tibility of the hlms can be increased within a range 
of relatively low No values that is high enough 
to pin domains in a head-to-head (or tail-to-tail) 
conhguration—this also requires sufficiently thick 
depletion layer thickness—but yet low enough to 
allow movement of domain walls with applied elec¬ 
tric helds. In our hlms, this range of No and w 
coincides with the formation of a wavy interface 
between head-to-head domains. Our results sug¬ 
gest that space charge-induced domains may engi¬ 
neer small hysteresis, large susceptibility response 
in ferroelectric materials for device applications. 
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